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ABSTRACT 

Context. The NO cycle takes place in the deepest layer of a H-burning core or shell, when the temperature exceeds T ^ 30 • 10* K. 
The O depletion observed in some globular cluster giant stars, always associated with a Na enhancement, may be due to either a deep 
mixing during the RGB (red giant branch) phase of the star or to the pollution of the primordial gas by an early population of massive 
AGB (asymptotic giant branch) stars, whose chemical composition was modified by the hot bottom burning. In both cases, the NO 
cycle is responsible for the O depletion. 

Aims. The activation of this cycle depends on the rate of the '^N(p,7)'*0 reaction. A precise evaluation of this reaction rate at 
temperatures as low as experienced in H-burning zones in stellar interiors is mandatory to understand the observed O abundances. 
Methods. We present a new measurement of the '^N(p,7)'*0 reaction performed at LUNA covering for the first time the center of 
mass energy range 70-370 keV, which corresponds to stellar temperatures between 65 -10' K and 780 -10* K. This range includes the 
'^N(p,7)'*0 Gamow-peak energy of explosive H-buming taking place in the external layer of a nova and the one of the hot bottom 
burning (HBB) nucleosynthesis occurring in massive AGB stars. 

Results. With the present data, we are also able to confirm the result of the previous R-matrix extrapolation. In particular, in the 
temperature range of astrophysical interest, the new rate is about a factor of 2 smaller than reported in the widely adopted compilation 
of reaction rates (NACRE or CF88) and the uncertainty is now reduced down to the 10% level. 

Key words, physical data and processes: nuclear reactions, abundances 

1 . Introduction stars the energy supply is provided by the pp-chairQ, the CNO 

cycle is the principal nuclear process in the core of high mass 
Hydrogen burning in stars proceeds through two diflFerent sets rnain sequence star s (M > 1.2 M p) as well as in the H-burning 
of nuclear reactions: the proton proton (pp) chain and the carbon shell of giant stars |Iben(1967)| Furthermore, a hot CNO cycle 
nitrogen oxygen (CNO) cycle. While in low mass main sequence may occur at the surface of H-accreting compact o bjects, hke 

white dwarfs or neutron stars |Jose & Hernanz(1998)| 
The set of nuclear reactions involved in the CNO cycle is 
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Fig. 1. The CNO bi-cycle. 



and the '^N(p,y)^^0, respectively: the ratio of the rates pro- 
vides the linkage between the CN and the NO cycles. The 
CN cycle becomes fully active when the temperature attains 
Tg > 0.016-0.02(3, while the NO cycle requires higher tem- 
peratures (Tg > 0.030-0.03^. In case of an active NO cycle, 
this process determines the abundances of all the stable oxygen 
isotopes ('^O, '^O, "*0). For this reason, a precise evaluation 
of the '^N(p,7)'*0 reaction rate is needed to address several 
astrophysical problems, like deep mixing scenarios in red gi- 
ant stars ?, see e.g.]sweigartl979, langerl986, charb98, kraft93, 
boo95, wass95, denis2003, pallO, hot bottom burning nucle- 
osynthesis in massive AGB stars |Renzini & Voli(1981) or the H- 
burning nucleosynthesis in nova-like events Iliadis et al.(2002)] 
IJose et al.(2007)l 

At low energies the cross section cr{E) of the '^N(p,y)"'0 
reaction (Q-value - 12. 127 MeV) is typically expressed in terms 
of the astrophysical S-factor 5 (£■) defined for this reaction as: 



S{E) = cr(£)£exp(212.85/ V£) 



(1) 



where E is the center of mass energy in keV. 

In hydrostatic H-burning, the Gamow peak energy of this 
reaction ranges between 30 and 100 keV. Larger values, up to 
300 keV, may be attained during explosive burning. In this en- 
ergy range, the astrophysical S-factor is influenced by two res- 
onances at £■ = 312 and 964 keVQ related to excited states in 
'''O at Ej, = 12440 and 130 90 keV, respectively The reaction 
rates reported in the NACR E [Angulo et al.(1999)| and the CF88 
ICaughlan & Fowler(1988) compilations are based on the direct 
measurement presented by Rolfs and Rodney (1974). However, 
more recent R-matrix studies Mukham edzhanov et al.(2008)[ 
|Barker(2008)| which also take into account a previous ANC 
measurement [Mukhamedzhanov et al.(2008)| suggested a sub- 
stantial reduction of the S(0) (i.e. the astrophysical factor at 
E = 0). This result is in agreement with older direct measure- 
ments|Hebbard(1 960)||Brochard et al.(1973)| 

This discrepancy prompted an in-depth study of the reaction 
at LUNA (Laboratory for Underground Nuclear Astrophysics). 
The LUNA facility has been designed to study nuclear re- 
actions of astrophysical interest at the same energies of 
the stellar interiors, by taking advantage of the ultra-low 
background ' Bemmerer et al.(2005)| |Caciolli et al.(2009)| of the 



^ Tg = T{K)/10\ 

' The activation temperatures of both the CN and the NO cycles 
depend on the actual amount of C, N and O nuclei and, therefore, on 
the stellar metallicity. 

* In the center of mass reference. Beam energies are given in the 
center of mass reference unless otherwise stated 



INFN-Gran Sasso underground laboratory (a detailed descrip- 
tion of LUNA and its experimental study of the pp chain 
and CNO cycle may be found in the following reviews: 
[Costantini et al.(20 09) Broggini et al.(2010)| l. First of all, a re- 
analysis of data taken with nitrogen gas target of natural isotopic 
composition (0.4% '^N) atE- 90-230 keV has been performed 
[Bemmerer et al.(2009)| Then, a new measurement has been car- 
ried out at LUNA and Notre Dame LeBlanc et al.(2010)[ HPGe 
detectors and enriched TiN solid targets have been used to cover 
a wide energy range, namely: E - 120-1800 keV. Although the 
minimum energy is still too high to study most of the stellar H- 
burning environments, thanks to the excellent accuracy (7%) and 
the wide energy range, this new experiment provided a dataset 
suitable for an R-matrix extrapolation toward lower energies. 

In this paper, we present a third experiment performed at 
LUNA, designed to explore lower energies. The use a BGO de- 
tector, having a higher y-detection efficiency compared to the 
HPGe detectors, allowed us to easily cover the 312 keV reso- 
nance region and to extend the direct measurements down to 70 
keV. The aim of this further effort is twofold. First of all, the new 
data set covers the Gamow peak corresponding to the explosive 
burning in Novae as well as hot bottom burning in massive AGB 
stars. Furthermore, it provides an independent test of the low en- 
ergy R-matrix extrapolation. 

In the next section we illustrate the experiment, the data anal- 
ysis and the results. In particular, a comparison of the present, in- 
dependent measurement with the low energy predictions of the 
R-matrix analysis |LeB lane et al.(2010)l leads to the conclusion 
that the '^N(p,7)"'0 reaction rate is now known within a 10% 
confidence interval. A summary of the astrophysical studies re- 
quiring an accurate evaluation of the '^N(p,y)"'0 reaction rate 
follows. 



2. The new underground experiment 

The target and the y-ray detection set-up are those used in 
previous measurements and have been already extensively de- 
scribed elsewhere (for instance see I Limata et al.(201(j)] ). The 
proton beam (30-150 //A) reaches the water cooled target af- 
ter passing a 5 mm diameter collimator and aim long copper 
tube, which is cooled to liquid nitrogen temperatures and works 
as a cold trap in order to prevent impurities scattered by the beam 
from depositing on the target surface. The pressure in the target 
chamber is 5 ■ 10"^ mbar and no carbon deposition on the target 
is observed after the irradiation. This is checked by performing 
scans of the profile of the '"'N(p,'x)'^0 resonance at Ep = 278 
keV. The target chamber works as a Faraday cup and provides 
the integral of the charge deposited, hence the average beam in- 
tensity, with an overall uncertainty of 2% (a -300 V high volt- 
age is applied to the cold trap to suppress the secondary electron 
emission). 

The target is surrounded by a 47r-BGO summing crys- 
tal (28 cm long, 20 cm diameter, and 6 cm coaxial hole, 
ICasella et al. (2002)1 ). The 47r-BGO is essential in order to in- 
crease the y-detection efficiency, which is calculated with a sim- 



ulation based on GEANT4 |AgostinelIi et al.(2003) and carefully 
checked with radioactive sources and with the y-ray produced 
by the proton induced reaction "B(p,y)'^C at the E = 149 keV 
resonance. The simulation needs experimental inputs, such as 
the decay scheme and the angular distribution of the emitted 
y-radiation. The decay branching ratios for transitions to the 
excited state of '^O have been measured by Rolfs and Rodney 
(1974), Bemmerer et al. (2009) and LeBlanc et al. (2010). The 
angular distribution has been found to be isotropic in a previous 
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Fig. 2. Schematic representation of the last portion of the beam- 
line with the detection set-up. 



LUNA work |LeBlanc et al.(2010)| By considering all the contri- 
butions described above in the simulation code, the total uncer- 
tainty on the efficiency is 3%. 

The TiN forming the target material, enriched in '""N, is 
deposited on a tantalum backing with the reactive sputter- 
ing technique |Rigato et al.(2001)| The target thickness is 100 
nm, as verified through secondary neutral mass spectroscopy 
|Vad et al. (2009)1 (the uncertainty on this measurement is in- 
cluded in the contribution to the target analysis in Table|2]l, corre- 
sponding to 15 keV energy loss at £ = 259 keV. The stoichiome- 
try Ti/N, which ranges from 0.97 to 1.18 according to the target, 
is measured for each target with the high Z elastic recoil de- 
tection (ERD) technique Bergmaier et al.(1998)| Isotopic abun- 
dances between 96% and 99%, according to the target, are de- 
duced from the observed height of the plateau in the yield of the 
'"^N(p,y)'^0 resonance at £ = 259 keV, and from the ERD data. 
The results from these two methods agree within 2%. Finally, 
the target deterioration, caused by the impinging high-intensity 
proton beam has been studied by using the 430 keV resonance 
of '^N(p,gy)^^C |Marta et al.(2010)| The targets have been an- 
alyzed by looking at the shape of the plateau in the yield dis- 
tribution for the 430 keV resonance. The surface irradiated by 
the LUNA beam and the area outside the LUNA beam-spot are 
investigated, so that appropriate corrections for the target deteri- 
oration during measurements are derived. 

The laboratory background in the '^N(p,y)"'0 region of in- 
terest is about 6 counts/day. The beam induced background in the 
same region, produced by the "B(p,')')'^C reaction, is monitored 
by means of the peak produced by this reaction at 16.1 MeV (see 
Fig.[3]as an example of the acquired spectra). The counts in this 
peak are usually more than in the 1 1 MeV peak produced by the 
same reaction, which Ues within the '^N(p,y)'^0 region of inter- 
est of the spectrum [Bemmerer et al.(2009)| We have rejected all 
measurements where the 16 MeV peak contained more than 3% 
of the counts in our region of interest. 

The target profiles can be integrated with the cross section in 
order to calculate the expected yield as: 






S (Ep) ■ Ep exp 



212.85 



p / 



15n 

'/BGO ■ nt^rgei{x)—dx (2) 



where Ep is the energy in the laboratory system expressed in keV 
and it depends on the beam position x along the target thick- 
ness, //BGO is the efliciency, and ntarget(^)-f|^ is the number of 
'^N nuclides in the x position in the target. By comparing the 
experimental yield Fgxp with the calculated one, it is possible to 
determine the S-f actor as follows: 




E^[MeV] 

Fig. 3. Spectra taken at fip = 330 keV (red), Ep = 80 keV (j 
and laboratory background spectrum (blue), normalized 
same time. 
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Table 1. The absolute S-factor data and their statistical uncer- 
tainties from the present work. The systematic error is 11%. 
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where the effective energy is calculated according to the follow- 
ing definition Lemut(2008): 



leff 



£- S(E)E exp (2>2|l) . n,,,,,,(x) ■ '-^ ■ E ■ dx 
£"'" S(E)-E exp (2i^) ■ n,,,get(^) ■ '-^dx 



(4) 



S (£'eff)exp - — ■ S (Eeff)i] 



(3) 



In Eq. (|4|i the theoretical S-factor is used. Four diff'erent the- 
oretical S-factors are considered in Eq. ^ and Eq. (HJl: the 
one reported in LeBlanc et al.(2010)l and the one reported in 
[Mukhamedzhanov et al.(201 1) a constant S-factor and a value 
obtained from a recursive analysis process. In all cases, the same 
results are obtained within 1 % discrepancies which is included 
in the error on the effective energy. 

As reported in Table [T] the '^N(p,y)'^0 astrophysical S- 
factor is obtained for the center of mass energy range 70-370 
KeV. The statistical uncertainty is always limited within a few 
percent, reaching a maximum value of 10% atE- 72.8 keV All 
sources of systematic uncertainties are given in Table|2]and sum 
to a total systematic uncertainty of 10%. 
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Table 2. S-factor systematic uncertainties. 



Source 


Estimated 


description 


uncertainty 


Target analysis 


7.5% 


stopping power 


4.0% 


'^N isotopic ratio 


2.0% 


Ti/N stoichiometry 


2.0% 


Beam intensity 


2.0% 


Effective energy 


3.0% 


y-ray detection efficiency 


3.0% 


"B(p,7)'2c background 


3.0% 


Totai systematic uncertainty 


10.0% 



Table 3. Summary of 3 12 keV resonance cross sections in com- 
parison to previous results (see text for details and references). 
The uncertainty reported by Hebbard (1960) has been obtained 
by assuming it to be 10% as reported by |Barker(2008)| 



A comparison of the derived astrophysical S-factor to the re- 
sults of previous experiments is shown if Fig.|4] We confirm the 
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Fig. 4. The S-factor as a function of energy. Present data (black 
dots) are compared to the results of previous experiments. 



previous finding concerning the need of a substantial reduction 
of the S(0) value. The present result is significantly lower than 
the resonant cross section from Rolfs and Rodney (1974), i.e. 
the data set adopted in NACRE and CF88, and, by considering 
the systematic uncertainty, in good agreement with our previ- 
ous HPGe measurement LeBla nc et al.(2010)| In particular, ac- 
cording to the present absolute analysis, the cross section on top 
of the £■ = 312 keV resonance is cr(312 keV) = 6.0 ± 0.6 //b, 
where the quoted error includes the 10 % systematic uncertainty. 
In Table [3] we compare this result to the values of previous 
measurements. The weighted average of 3 measurement^ leads 
to a recommended value of (t(312 keV) - 6.5 ± 0.3 fih. The 
shape of the R-Matrix fit has been also compared to the present 
data as shown in Fig. |5] Only for this comparison the present 
data have been corrected for the electron screening in the adi- 
abatic approximation Asse nbaum et al.(1987)| (at most 10% at 
70 keV) and they have been rescaled to the calculated average 
value. This rescaling is still between inside the systematic un- 
certainties of the present absolute data. They show an excellent 
agreement with the energy dependence of the LUNA R-Matrix 
fit |LeBlanc et al.(2010)| 



^ The result obtained by |Brochard et al.(1973)| has been excluded, be- 
cause no uncertainty was reported. 
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Fig. 5. The present data rescaled to the new recommended value 
of cr(3 12) keV (see text to details) and compared to the R-matrix 
predictions [LeBlanc et al.(201"0)| 



Finally, a new R-matrix analysis has been recently published 
by Mukhamedzhanov, La Cognata and Kroha (2011). By vary- 
ing the fitting method, these authors obtain S (0) values ranging 
between 33.1 and 40.1 keVb, which is in excellent agreement 
with the value reported by LeBlanc et al. (5(0) = 39.6 + 2.6 
keVb). 

For practical purposes, the nuclear reaction rate 
can be approximated by the following fitting formula 
[LeBlanc et alX20TT)| 



Na<o-v> ^ ai 10V"5 exp[a27'"^ - (T/oif] (5) 

[1 +a4T + a5T^]+a6lO^T--2 
expiai/T) + aglO^T-J exp(fl9/r), 

where the best fit parameters are reported in Table|4] 



Table 4. Best fit parameters for the '^N(p,y)'^0 reaction rate 
given in [LeBlanc et al.(2010)[ 



ai = 0.523 


a4 = 6.339 


ay = -2.913 


a2 = -15.240 


as =-2.164 


as = 3.048 


aj = 0.866 


ag = 0.738 


a, = -9.884 



3. Summary and Conclusions 

In this paper we have discussed the experimental eff'orts done to 
improve our knowledge of the '^N(p,y)'^0 reaction rate in the 
temperature range experienced by any H-burning zone in stellar 
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interiors. Such an important reaction is located at the branch- 
ing point between the CN and NO cycles. The branching ratio, 
as a function of the temperature, is shown in Fig. |6l where the 
solid line has been obtained by means of the widely adopted re- 
action rate given by NACRE, while the dashed line represents 
the revised scenario as derived from the latest R-matrix study 
(see section|2]l. In both cases, the rate suggested by NACRE has 
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Fig. 7. The '^O equilibrium abundance (mass fraction) as given 
by: Zi6 = Zisll ^'^'^''"^'' . The solid and the dashed lines repre- 

sents the old (NACRE) and the new (revised) prescriptions for 
the '^N(p,'x)'*0 reaction rate. A solar '^N mass fraction (X15), 
as derived by Asplund et al. 2010, has been used. In both cases, 
NACRE prescriptions for the < crv >i6o+p have been adopted. 



Fig. 6. The CN-NO branching ratio, as a function of the temper- 
ature, under different assumptions for the '^N(p,y)'^0 reaction 
rate: NACRE (solid line) and the revised rate (dashed line). 

been used for the competitive '^N(p,a)'^C reactiorQ A look at 
the solid line shows that in the whole range of temperatures ex- 
perienced by the core and the shell-H burning, the a-channel 
is between 1000 to 2000 times more efficient than the y chan- 
nel: just 1 to 2 protons out of every 2000 are consumed by the 
NO cycle. When the updated rate for the '^N(p,7)'^0 is adopted, 
such a ratio becomes about a factor of 2 larger. Although such 
a variation has negligible consequences on the overall nuclear 
energy production, a change in the rate of the '^N(p,'y)^^0 af- 
fects the equilibrium abundances of the stable oxygen isotopes 
within the H burning zone. As an example, the equilibrium abun- 
dance of "'O is reported as a function of the temperature in Fig. 
I2I Also in this case, the solid and the dashed lines represent the 
values obtained by adopting the NACRE and the revised rate of 
the '^N(p,y)'^0 reaction, respectively. 

Let us point out that the most important improvement re- 
sulting from the present analysis of the CN-NO branching con- 
cerns the significant reduction of the nuclear physics uncertain- 
ties, other that the change of the reaction rate with respect to the 
values reported by CF88 or NACRE. For stellar models and nu- 
cleosynthesis calculations implying H-burning whose Gamow 
peak energy is larger than the minimum value attained by the 
LUNA BGO experiment, namely £0 > 70 KeV, which corre- 
sponds to a temperature T > 65 ■ 10'' K, a true experimental 
eiTor (smaller than 10%) is now available for this important re- 
action rate. Note that only in a very few cases the reaction rate 

* This reaction has been recently studied with the THM method 
|Cognata et al.(2 009) The authors do not report a reaction rate but only 
the 5(0) value. Scaling the previous NACRE results on that value the 
following considerations do not change so we still adopt the NACRE 
results in the present work 



has been measured down to the stellar Gamow peak energy (see, 
e.g., |Bonetti et al.("l999)) . In addition, basing on the good agree- 
ment found between the new LUNA measurements and the re- 
vised R-matrix fit (see previous section), we are confident that 
the quoted small uncertainty may be assumed also in the extrap- 
olated region. 

Among the many astrophysical applications of the present 
analysis, we recall the explosive H-burning in Novae, which oc- 
curs at temperature larger than 10** K and, therefore, well above 
the achieved experimental limit. A recent study by Iliadis et al. 
(2002), investigates the dependence of the nova nucleosynthesis 
calculations on the various nuclear physics inputs. They found 
that a reduction of a factor of two of the '^N(p,7)"'0 reaction 
rate would imply a 30% reduction of the final oxygen abundance. 
Also the inner region of the convective envelope of massive 
AGB stars attains quite high temperature, up to Tg = 0.08-0.09 
'Renzini&Voli(1981)' 'M. Forestini & C. Charbonnel(I997) 
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D'antona & Mazzitelli(1996^ 

Lattanzio et al.(2000) The resulting H burning, the so called 
hot bottom burning, coupled to the convective mixing, gives rise 
to a very promising nucleosynthesis scenario, where all the C, 
N and O isotopes are substantially affected. If the temperature is 
large enough (80 x 10^ K), the Ne-Na and the Mg-Al cycles are 
also activated. In this context, it has been recently claimed that 
massive AGB stars played a fundamental role during the early 
evolution of globular clusters [Ventura et al.(200T)| According 
to this self-enrichment scenario, in between 50 to 100 Myr 
after the cluster formation, the first generation of intermediate 
mass stars (5-7 Mq) reached the AGB. Then, during this 
evolutionary phase, they underwent a substantial modification 
of the envelope composition, as a consequence of the HBB 
and several dredge up episodes. Due to the huge AGB mass 
loss, fresh gas enriched in He, C, N and Na, but O depleted, 
refilled the space occupied by the young Globular Cluster. If 
the star formation process was still active at that epoch, some 



A. CacioUi et al.: Revision of the "NCp,-/)'*© reaction rate and oxygen abundance in H-burning zones 



of the stars we observed today should show the imprint of such 
a delayed chemical pollution by massive AGB. In particular, 
the O-Na anti-correlation, as observed in Giant, sub-Giant and 
turn-off stars of several globular clusters ?, e.g. and reference 
therein]kraftl997, Carretta2009, may be the consequence of 
this nucleosynthesis process. Such a conclusion follows from 
the evidence that the temperature required for the activation 
of the NO cycle is similar to that required for the activation 
of the Ne-Na cycle. Thus, when O is depleted at the bottom 
of the convective envelope, Na should be enhanced. For this 
reason, a precise determination of the '^N(p,7')"'0 is one of the 
prerequisites to obtain a robust prediction of the O abundance 
and, in turn, to check the proposed self-pollution scenario for 
the observed O-Na anti-correlation. 

The R-matrix studies also allow to extrapolate the precise 
experimental measurements of the '^N(p,7)'^0 reaction rate 
down to the temperature range experienced by the H-burning 
taking place in main sequence, RGB and less-massive AGB 
stars. Also in these cases the uncertainty has been significantly 
reduced. Such an occurrence may be immediately translated in 
more robust astrophysical predictions. 
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